Electron Teleportation via Majorana Bound States in a Mesoscopic Superconductor 
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Zero-energy Majorana bound states in superconductors have been proposed to be potential build- 
ing blocks of a topological quantum computer, because quantum information can be encoded in 
the fermion occupation of a pair of spatially separated Majorana bound states. However, despite 
intensive theoretical efforts, non-local signatures of Majorana bound states have not been found 
in charge transport. In this work, we predict a striking non-local phase-coherent electron transfer 
process by virtue of tunneling in and out of a pair of Majorana bound states. This teleportation 
phenomena only exists in a mesoscopic superconductor because of an all-important but previously 
over-looked charging energy. We propose an experimental setup to detect this phenomena in a 
superconductor/quantum spin Hall insulator/magnetic insulator hybrid system. 

PACS numbers: 71.10.Pm, 73.23-b, 74.50+r, 74.90.+n 



Majorana bound states are localized zero-energy exci- 
tations of a superconductor QJ An isolated Majorana 
bound state is an equal superposition of electron and hole 
excitations and therefore not a fermionic state. Instead 
two spatially seperated Majorana bound states together 
make one zero-energy fermion level [TJ |3] which can be 
either occupied or empty. This defines a two-level sys- 
tem which can store quantum information non-locally. as 
needed to realize topological quantum computation[31[5]. 
While several schemes have been recently proposed to de- 
tect the existence of individual Majorana bound states 0- 
112] , experimental signatures of the non-local fermion oc- 
cupation of these states remain to be found. 

In this work, we predict a non-local electron transfer 
process due to Majorana bound states in a mesosopic 
superconductor: an electron which is injected into one 
Majorana bound state can go out from another one far 
apart maintaining phase coherence. Strikingly, the trans- 
mission phase shift is independent of the distance "trav- 
elled". In such a sense, we call this phenomena "elec- 
tron teleportation". It occurs because of the non-local 
fermion occupation of Majorana bound states and the 
finite charging energy of a mesoscopic superconductor. 
The all-important role of charging energy in the study 
of Majorana fermions has not been recognized before. 
We propose a realistic scheme to detect the teleportation 
phenomena in a superconductor/quantum spin Hall insu- 
lator/magnetic insulator hybrid system, which have been 
recently shown to host Majorana bound states [HAITI]. 

In a macroscopic s-wave superconductor, charge e ex- 
citations have a pairing energy gap, whereas charge 2e 
excitations cost zero energy. Therefore the ground state 
manifold consists of states with an even number of elec- 
trons only, as shown in Fig. la. The BCS wave-function of 
the ground state with a definite overall superconducting 
phase <j) € [0, 2n], is a linear superposition of states with 
2N electrons. Now consider that a pair of zero-energy 
Majorana bound states are present at positions R\ and 
i?2 in the superconductor, and all other quasi-particle ex- 
citations have a finite gap greater than an energy scale A. 




FIG. 1: Energy spectrum of a superconductor as a function 
of total number of electrons. States with an even and odd 
number of electrons are marked in black and red respectively, 
(a) and (b) correspond to superconductors without and with 
a pair of zero-energy Majorana bound states. Figures on the 
left and right correspond to superconductors without and with 
charging energy. 

We shall show later how this situation can be realized in 
a device consisting of an s-wave superconductor and the 
recently discovered quantum spin Hall (QSH) insulator 
HgTe quantum well[T3[TS]. The two Majorana operators 
71 and 72 are defined by: 

71,2 = J dxe-Wfaixtfix) + e^ 2 C ia {x)c{x). (1) 

Here ^1,2(2^) are bound state wave-functions centered at 
Ri.2- We assume that the distance between the two 
Majorana bound states is much larger than the coher- 
ence length — a necessary condition for the notion of non- 
locality to be meaningful. A single fermionic operator 
can then be defined d) = (71 + i72)/2, which accom- 
modates an extra fermion excitation without any energy 
cost. Now the ground states of the superconductor have 
two sectors |e) and |o) defined by d\e) = 0, |o) = d)\e), 
which have an even and odd number of electrons respec- 
tively as shown in Fig. lb, 
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where a n is real and slowly varying at large n. Eq.([2]) 
says that the fermion occupation of the d level (empty 
or occupied) (fid — (17172 + l)/2, is /raec? by the total 
number of electrons in the superconductor mod 2: 



«7i72 = (-I)™- 



(3) 



Eq.([3]) imposes constraint on the Hilbert space. Equiv- 
alcntly, ^ implements a gauge transformation jj — > 
—7,-, (f> — > 4> + 27T, which is a gauge symmetry in the def- 
inition of Majorana operators Eq.Q. 

If the superconductor under consideration is of meso- 
scopic size and connected to ground by a capacitor, the 
energy spectrum has an additional term due to the finite 
charging energy: 



U c (n) 



Q ) 2 /2C, n = 0,±l,±2, 



(4) 



where Qo is the gate charge determined by the gate volt- 
age V g across the capacitor. As a result, states with 
different n are no longer degenerate. In this work we 
will consider the regime U = e 2 /C < A, which can 
always be satisfied by increasing the size of the super- 
conductor. The low-energy spectrum (E < U) then de- 
pends crucially on whether Majorana bound states are 
absent or present. In the former, only states with an even 
number of electrons are important at low energy, which 
leads to an even-odd effect in tunneling experiments on 
mesoscopic superconductorsfTT] . In contrast, if Majorana 
bound states are present, both even- and odd-states ap- 
pear on equal footing in the low-energy spectrum. In this 
case, when Qo/e is adjusted to half-integers, an energy- 
level degeneracy can be achieved between two states that 
differ by charge e, instead of 2e as in a usual supercon- 
ductor with Coulomb blockade. This two-level system 
shown in Fig.l is the main subject of our study. 

We now weakly couple the two Majorana bound states 
to separate normal metal leads. This can be realized 
in an s-wave superconductor/quantum spin Hall insula- 
tor/magnetic insulator hybrid system. The QSH insula- 
tor used here is a new phase of two-dimensional insulators 
which have robust edge states [IS]. It has been experi- 
mentally realized in HgTe quantum wells [TS]. The device 
geometry is shown in Fig. 2: an s-wave superconductor 
and a magnet are deposited on top of the QSH insulator. 
Both proximity effect with the superconductor and the 
Zeeman field due to the magnet open up a finite quasi- 
particle gap for the QSH edge states. However, two zero- 
energy Majorana bound states 71,2 exist at the intersec- 
tion of the superconductor-magnet interface with the top 
and bottom edge respectively [Mj, conceptually similar to 
the states localized at the ends of a one-dimensional spin- 
less p-wave superconductor [TU]. The edge states in the 
uncovered part of the quantum spin Hall insulator are 
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FIG. 2: the device used to study electron tunneling from leads 
into two Majorana bound states, consisting of an s-wave su- 
perconductor (SC) and a magnetic insulator (M) junction de- 
posited on top of a quantum spin Hall insulator (QSH). By 
tuning the gate voltage across the capacitor, the supercon- 
ductor is close to a charge degeneracy point between a total 
number of electron no and (no + 1). This two-level system cor- 
responds to a resonant level empty or occupied. At a small 
bias, electron tunneling through two Majorana bound states 
is equivalent to phase-coherent tunneling through a single res- 
onant level. 



naturally used as leads to connect 71 and 72 to source 
and drain. 

To describe electron tunneling between the lead and 
the superconductor, we write the electron operator in 
terms of quasi-particle operators of the superconductor 



c(x) = e-^' 2 [^{x) ll+ i 2 {x) l2 + 



(5) 



Since we will only consider small bias voltage V < U < 
A, only zero-energy Majorana operators are important 
and contributions from other quasi-particle states can be 
neglected in ([5]). We now write down the Hamiltonian 
for the system in Fig. 2: H = Hi + U c + Ht, where Hl = 
Sfc j=i 2 e 3 : (^) c j k c j,k describes the two leads, U c is the 
charging energy defined in Q. The effective tunneling 
Hamiltonian at low energy is given by substituting 
(|5| into the bare tunneling term ^eje: 



Ht = 



E 

J=l,2 



0/2 



(6) 



Here Cj annihilates an electron in lead j and Xj oc £,j(Rj) 
is the tunneling matrix element. As we emphasized ear- 
lier, £,i(x) and £2(2:) have essentially zero wavefunction 
overlap so that no coupling between €1(02) and 72(71) 
exists. The operator e ±J ^/ 2 in Ht increases/decrases the 
total charge of the superconductor by one charge unit 
[n, e ±l ^/ 2 ] = ±e ±41 ^/ 2 , and the Majorana operator 7^2 
changes the parity of electron number in the supercon- 
ductor. The "naive" Hilbert space of H is simply the 
direct product of electron number eigenstate \n) and the 
two states of d- level (|e) and |o)), but it is redundant. 
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Instead the physical Hilbert space only consists of those 
states \n = 27V; e) and \n ~ 27V + 1; o) that satisfies the 
gauge constraint 

When the source is biased at a small voltage V, cur- 
rent flows to drain by electron tunneling in and out of the 
superconductor. Since charging energy U c favors states 
with a fixed number of charge in the superconductor, 
only two charge states |no) and \ng + 1) give dominant 
contribution to the current for V < U, which is simi- 
lar to tunneling through a quantum dot in the Coulomb 
blockade regime. To a good approximation, we can then 
truncate the Hilbert space only keeping these two states, 
which we label by s z — ±1. H then becomes 



H 



.7=1,2 



(7) 



Here 5 is the energy difference between |no) and \hq + 1) 
and depends on the gate voltage. The gauge symmetry 
^ then becomes ijij 2 s z — (—1)"°. The key to solve the 
tunneling problem Q is to combine Majorana and spin 
operators into a singe fermion operator /. 



lis + -> / + , 7i«~ 
l2 s+ -> i(-l) n °f + , l2 s~ 



f 

i (_l)«o+i / . (8) 



We have checked that this transformation preserves all 
commutation relations 



>7j s+ } = {7^",7iS^} = 0, 
{7is + ,72S~} - 7172S, =i(-l) no+1 , 



(9) 



where the gauge constraint is used in the last equation. 
Conceptually, it is not surprising that the transforma- 
tion Q works: after all, the two charge states |no) and 
|n + 1) differ by one electron. Using ([8]), we rewrite the 
Hamiltonian H in terms of the fermion operator /: 



H= H L + S(pf-^) + (X lC \f + hc) 
+(-l) n °(-i\ 2 clf + hc). 



(10) 



Eqs. ([6j [7j [8j [To]) are the main results of this work, and to 
the best of our knowledge they have not been reported 
before. ( fl0| ) says that electron tunneling in and out of two 
spatially separated Majorana bound states is equivalent 
to resonant tunneling through a single level, as shown 
schematically in Fig. 2. Since resonant tunneling is a co- 
herent process, we conclude that an incident electron at 
E < U tunnels into one Majorana bound state and comes 
out from its partner far apart still maintaining phase 
coherence. Strikingly, the magnitude and phase of the 
transmission amplitude — which we call t — is independent 
of their distance. In this sense, we call such a non-local 
electron transfer process "teleportation" . 
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FIG. 3: Left: an interferometer that probes the phase- 
coherent electron teleportation via two Majorana bound 
states. Right: schematic plot of zero bias differential conduc- 
tance as a function of magnetic flux at two successive charge 
degeneracy points (27V - 1, 27V) and (27V, 27V + 1). The h/2e 
shift in conductance peak signals the change of fermion num- 
ber parity in the superconductor. 



The phase coherence over a long distance shown here 
is in fact a direct consequence of Majorana bound states. 
Conceptually it can be best understood from electron's 
Green function G e ^ (x, t; y, 0) = (c(x, t)c^(y, 0)) e ( o ) de- 
fined in the even and odd ground state sector |e) and |o) 
respectively [20 . Using ([5]), we find 

G e >°(x,t^oo;y,0) = TtQ(v)ti{x) ~ 0(1) (11) 

is finite for x ~ Ri, y ~ R 2 - The long-time limit corre- 
sponds to the low-bias regime we are interested in. The 
fact that (11 1 is independent of — R 2 \ is most unusual, 



as first pointed out in Ref . |2"T] . 

We now show that, interestingly, the phase of trans- 
mission amplitude t depends on the gate charge Qo in a 
surprising way, and therefore it is sensitive to the fermion 
occupation (|e) vs. |o)) of the Majorana bound state pair. 
Consider tuning gate voltage to make Qo change by one 
charge unit. The number of electrons in the supercon- 
ductor ground state will correspondingly change by one. 
Although the excitation energy spectrum U c comes back 
to itself, we find 



t 



-t, when Qo Qo ± e, 



(12) 



i.e., the transmission phase shift changes by n. This be- 
havior is related to the change of fermion number parity 
in the ground state. The property ( 12 ) is evident from 
the (— l) n ° factor in H, which valid in the two-level ap- 
proximation. Using second-order perturbation theory in 
the weak tunneling limit, one can easily show that this 
(—1)™° factor comes from the ± sign in Eq.([lT]). In gen- 
eral, we can prove (12 1 using the following symmetry of 
the full Hamiltonian H and the gauge constraint Q 



C/ J ff(Q ,A 2 )^ 1 
C/j 7l72 (-l)«f/-i 



H(Q + e, -A 2 ), 
«7i72(-l)", U = 



72 e 



-10/2 



To detect the phase-coherence of the electron telepor- 
tation described above, we consider the interferometer 
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FIG. 4: Compared with Fig. 2, the superconductor here is 
grounded without charging energy. Charge transfer between 
the lead and the superconductor is conducted by local An- 
dreev reflection that transfers charge 2e. 



setup in Fig. 3: a point contact is introduced to partially 
scatter an incident electron at the top edge directly to 
the bottom edge, and partially transmit it to the super- 
conductor which subsequently comes out at the bottom 
edge. Interference between the two paths with a mag- 
netic flux <I> enclosed leads to a ^-dependent differen- 
tial conductance dl/dV, which is /i/e-periodic. If the 
direct scattering probability is large, interference visibil- 
ity is maximum at the charge degeneracy point when 
electron tunneling through Majorana bound states is on 
resonance. We schematically plot the <E> dependence of 
dl / dV (at zero bias and zero temperature) for two suc- 
cessive half-integer charges Qq in Fig. 3. The sign reversal 



of t discussed in ( 12 ) leads to a h/2e shift in the interfer- 
ence pattern. 

Discussion: It is instructive to compare our study 
of tunneling into Majorana bound states in the V < 
U regime with previous studies which do not include 
the charging energy U 6 8 . Instead of using a "float- 
ing" superconductor as in Fig. 2, these works consider a 
grounded superconductor in tunneling contact with two 
leads each having an independent bias voltage. Such a 
three-terminal device can be realized in a geometry shown 
in Fig. 4. In this setup, transferring two electrons to the 
superconducting condensate does not cost energy. There- 
fore, when the two Majorana bound states are sufficiently 
far apart, an incident electron from a lead can be An- 
dreev reflected as a hole to the same lead, but will never 
appear in the other lead. In other words, no electron 
teleportation happens there. Indeed, Bolech and Dern- 
ier has shownjf)] that: (a) the conductance of each lead 
is 2e 2 /h at zero bias and zero temperature, which is a 
sign of charge 2e transfer by Andreev reflection; (b) the 
two leads have independent currents without any cor- 
relation (no teleportation). The same results were also 
obtained using scattering approach within Bogoliubov-de 
Gennes formalism[5]. In the two-terminal device we stud- 
ied (Fig. 2), charge 2e transfer is suppressed by charging 
energy and the conductance is at most e 2 /h because of 
single electron tunneling. 

Finally, we discuss the feasibility of experiments using 



the quantum spin Hall insulator HgTe quantum well. A 
good candidate for the superconductor in our proposed 
setup is indium with T c = 3.4K, which is currently used 
as electrodes to contact HgTe [22]. The relevant parame- 
ters for these materials have been estimated [5] [H]. As- 
suming a proximity- induced gap A ~ O.lmeV, the pen- 
etration length of Majorana bound states is about 3^tm. 
So if the top and bottom Majorana bound states in Fig. 2 
are 30/^m apart, direct tunneling between them are neg- 
ligible. If charging energy of the superconductor can be 
optimized to be comparable to A, the resonant-tunneling 
model for electron teleportation described above is valid 
below IK. The level broadening T from coupling to 
leads sets the temperature scale for detecting the phase- 
coherence. 

In summary, we reveal a striking non-local electron 
transport phenomena through Majorana bound states in 
a finite-sized superconductor with charging energy. Most 
interestingly, the transmission phase shift detects the 
state of a qbit made of two spatially separated Majorana 
bound states. In a future work [23]. we will propose a gen- 
eralized scheme to electrically detect the internal states 
of multiple Majorana bound states, with an emphasis on 
implementing topological quantum computation. 
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